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We analyze the steady radial accretion of matter into a nonrotating black hole. Neglecting the self¬ 
gravity of the accreting matter, we consider a rather general class of static, spherically symmetric and 
asymptotically flat background spacetimes with a regular horizon. In addition to the Schwarzschild 
metric, this class contains certain deformation of it which could arise in alternative gravity theories or 
from solutions of the classical Einstein equations in the presence of external matter fields. Modeling 
the ambient matter surrounding the black hole by a relativistic perfect fluid, we reformulate the 
accretion problem as a dynamical system, and under rather general assumptions on the fluid equation 
of state, we determine the local and global qualitative behavior of its phase flow. Based on our 
analysis and generalizing previous work by Michel, we prove that for any given positive particle 
density number at infinity, there exists a unique radial, steady-state accretion flow which is regular 
at the horizon. We determine the physical parameters of the flow, including its accretion and 
compression rates, and discuss their dependency on the background metric. 

PACS numbers: 04.20.-q,04.70.-s, 98.62.Mw 


I. INTRODUCTION 

The accretion of gas into a black hole is an interesting physical process which has been studied extensively in the 
literature. As the gas falls into the black hole, effects due to compression or friction give rise to the emission of 
electromagnetic radiation whose spectrum can be observed and compared with theoretical models [I]. A distinctive 
feature of black hole accretion (as opposed to accretion into a star) is the presence of an event horizon acting as 
a one-way membrane through which the accreting gas disappears. This has several important implications. First, 
the horizon provides a natural inner boundary condition for the equations of motion describing the gas, and thus 
complications or uncertainties related to finding the correct boundary conditions at the star’s surface are avoided. 
Next, radiation emitted in the vicinity of the black hole undergo a strong gravitational lensing effect, resulting in an 
image showing the black hole’s shadow surrounded by a sharp light ring. In fact, millimeter-wave very-long baseline 
interferometric arrays such as the Event Horizon Telescope [2] are already able to resolve the region around Sagittarius 
A*, the supermassive black hole lying in the center of our galaxy, to scales smaller than its gravitational radius [3], 
and thus may lead to new tests of general relativity in its strong held regime, see for example Ref. [4]. 

In this article, we consider the simple case in which the black hole is nonrotating and the how is radial and steady- 
state. The problem of spherical accretion of a huid onto a gravity center has hrst been addressed by Bondi [51 
in the Newtonian limit. Generalizations to relativistic spherical steady-states on the Schwarzschild geometry and 
applications to the accretion of matter from the interstellar medium into a nonrotating black hole have been given by 
Michel [Bj and further discussed in Ref. (Tj . Given a huid equation of state and the temperature of the gas at inhnity, 
regularity of the how at the black hole horizon singles out a unique solution which turns out to be transonic, the sonic 
sphere describing the transition of the how’s radial velocity measured by static observers from sub- to supersonic 
lying outside the horizon. The stability of the Michel how in the subsonic region outside the sonic sphere has been 
established by Moncrief [7j. For more recent work on radial hows, see Ref. [8] for a perturbative study, Refs. EHm 
for the study of radial accretion onto cosmological black holes, and Ref. m for the accretion of dark matter into a 
supermassive black hole. For the self-gravitating case, see Refs. mm- 

The purpose of this work is to provide rigorous results regarding the existence and uniqueness of the Michel how, 
and to generalize it to equations of state more general than the simple polytrope and to background metrics more 
general than the Schwarzschild one. In fact, we consider a rather general class of static, spherically symmetric 
and asymptotically hat spacetimes which have a regular Killing horizon and whose metric coefficients obey certain 
monotonicity conditions. As a consequence, our results are also applicable to nonrotating black holes predicted by 
alternative theories of gravity or to distorted black holes in Einstein gravity in the presence of external helds, like 
dark matter for example. 

We study the accretion problem by casting it into a Hamiltonian dynamical system on the phase space (r, n) 
parametrized by the areal radius coordinate r and the particle density n of the huid. Although the Hamiltonian 
describing this system is hctitious and does not have the standard “kinetic plus potential energy” form, nevertheless 
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the system can be analyzed using standard tools from the theory of dynamical systems. Based on our assumptions on 
the fluid equation of state and the metric fields we give a complete description for the qualitative behavior of the phase 
flow. In particular, we show that there exists a unique hyperbolic critical point of this system which corresponds to a 
saddle point of the Hamiltonian. The global behavior of the stable and unstable one-dimensional manifolds through 
this critical point are discussed, and we prove that the unstable manifold describes the unique flow which extends 
from the horizon to infinity and has a finite, positive particle density at infinity. This solution describes the Michel 
flow and the critical point through which it passes corresponds to the sonic sphere. Conversely, we prove that for a 
given positive particle density at infinity there exists a unique critical point of the system whose unstable manifold 
matches it. Finally, assuming a polytropic equation of state with adiabatic index 7 lying in the range 1 < 7 < 5/3 
and focusing on the case where the fluid’s sound speed at infinity is much smaller than the speed of light c, we 
compute the accretion rate as a function of the particle density at infinity and the compression rates of the fluid at 
the horizon and the sonic sphere. For a gas with l<7<5/3we find that the compression rate at the sonic sphere 
does not depend on the details of the background metric to leading order in Uqo/c, though a nontrivial dependency 
appears to next-to-leading order. However, the compression rate at the event horizon depends quadratically on the 
ratio (7 1 = r s /rjj , where r s = 2 GM/c 2 is the Schwarzschild radius of the black hole and ru its event horizon radius 
(tii = 1 in the Schwarzschild case). For the case where the adiabatic index is 7 = 5/3 both the compression rates at 
the sonic sphere and the horizon depend on the details of the metric to leading order in Uoo/c. The accretion rate, 
on the other hand, is independent of the details of the background metric and is given by the familiar Bondi formula 
which only depends on the total mass M of the black hole, the particle density at infinity and Uqo/c to leading order 
in Voo/c. 

The remaining part of this paper is structured as follows: in Sec. |H] we first define the class of static, spherically 
symmetric black hole metrics we are considering in this work. Then, we review the derivation for the equations of 
motion describing a steady-state, radial perfect fluid flow, specify our assumptions on the fluid equation of state and 
enunciate our main results. In Sec. |III| we reformulate the fluid equations as a (fictitious) two-dimensional Hamiltonian 
dynamical system on the phase space (r, n), where the Hamiltonian depends on the accretion rate /i. We perform a 
detailed analysis of the phase flow, starting with the standard local analysis, and show that under our assumptions 
and for large enough |/u| there exists a unique hyperbolic critical point ( r c ,n c ), with associated stable and unstable 
local one-dimensional manifolds. Next, we analyze the global behavior of the phase flow and prove that the unstable 
manifold extends from the horizon to the asymptotic region, assigning to each value of r a unique value for the particle 
density n. This solution describes the Michel flow, and the critical point it passes through corresponds, physically, 
to the sonic sphere. We also show in Sec. |III| that given a positive particle number 77.00 > 0 at infinity, there exists 
a unique value for \/i\ and a unique critical point (r c ,n c ) such that the corresponding Michel flow solution matches 
the given value of 77 ^ as r -> 00 . We also make brief comments on other solutions which could be useful when 
describing stellar winds or subsonic flows accreted by relativistic stars. Next, in Sec.[TV]we discuss specific examples, 
including the accretion of dust and of a stiff fluid, provide numerical plots for the phase flow of polytropic fluids on a 
Schwarzschild and other backgrounds, and analyze the parameters characterizing the accretion flow under the usual 
assumption that the speed of sound at infinity is much smaller than the speed of light. For the case of a polytropic 
equation of state, we compute the accretion rate and the compression rates of the fluid at the sonic sphere and the 
event horizon for small Uqo/c, and analyze their dependency on the background metric. Conclusions are drawn in 
Sec. [VJ Some physical justification for our assumptions on the metric fields are discussed in Appendix [Aj while the 
proofs of some technical results needed in our analysis can be found in Appendix [B] 


II. STEADY-STATE RADIAL ACCRETION FLOWS INTO A BLACK HOLE AND MAIN RESULTS 

In this section, we specify our assumptions on the background black hole metric and on the fluid equation of state, 
review the relevant equations describing a relativistic, steady, radial flow, and state our main results regarding the 
existence and main properties of such flows. 


A. Assumptions on the metric field 

The background metric on which we study the accretion in this work is modeled by a spherically symmetric metric 
of the form 


ds 2 = —a(r) 2 N(r)c 2 dt 2 + 


dr 2 
N(r) 


r 2 (d’d 2 + sin 2 fldip 2 


(1) 
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where a and N are smooth functions of the areal radius r which are independent of the time coordinate t and the 
spherical coordinates ($, ip). Accordingly, the spacetime described by the metric ([I]) is invariant with respect to the 
flow defined by the Killing vector field k = 1 and with respect to rotations on the two-sphere. The geometric 
interpretation of the functions a and N is the following: the Killing vector field k and the differential dr of the areal 
radius give rise to the two scalar quantities a := —g(k, k) = —g llv k ,l k v and g (dr, dr) = 5 /i!y (V M r)(V l/ r). In terms of 
these scalar fields, N = g(dr, dr) and a = y/cr/N. For the Schwarzschild metric, a = 1 and N(r) = cr(r) = 1 — r#/r, 
and the event horizon is located at r = r#, where rjj = 2 GM/c 2 with M the black hole’s mass and G Newton’s 
constant. 

In this work, we consider a more general class of static, spherically symmetric black holes which is characterized by 
the following requirements on the smooth functions a, N : (0, oo) —> R: 

(Ml) Asymptotic flatness with positive total mass: 

As r —> oo, a(r) 1, N(r) —> 1, and r 2 a'(r) —> 2m for some positive constant m > 0. 

(M2) Regular Killing horizon: 

There exists a radius rn > 0 such that a(rn) = 0, er'(rf/) > 0 and a(ru) > 0. 

(M3) Monotonicity conditions: 

The functions a , a and a' are strictly positive on the interval r > th- 

(M4) Proximity to Schwarzschild condition: 

The following inequality holds for all r > r#: 


3 < 


(rVy 

ra' 


m 

< — 
r 



( 2 ) 


Before we proceed, let us briefly discuss the meaning of these conditions. For more details, definition and exam¬ 
ples we refer the reader to Appendix [A] Condition (Ml) assures that the spacetime described by the metric ([T]) is 
asymptotically flat in the sense that it approaches the Minkowski metric as r —> oo, and that its total Komar mass 
is M = mc 2 /G > 0. Next, condition (M2) implies that the metric possesses a Killing horizon at r = r#, where the 
Killing vector field k is null. The associated surface gravity is 


c 2 cr'(r H ) 
2 a(r H ) 


> 0 , 


(3) 


implying that the horizon is non-degenerated. Next, condition (M3) can be motivated by requiring the effective 
stress-energy tensor obtained by plugging the metric in Eq. |T]) into Einstein’s equations, to satisfy the weak and the 
strong energy conditions. Finally, for the meaning of the last condition we note that for the Schwarzschild metric, 
<j(t) = 1 — rn/r such that r 2 a' = r# is constant, implying that the inequality in Eq. (|2|) is trivially satisfied in this 
case. Therefore, the condition (M4) can be interpreted as a restriction on the size of deformation of the function a 
relative to its Schwarzschild value. An explicit example for such a deformation is given in Appendix [A] 


B. Radial perfect fluid flows 

After having specified our assumptions on the background geometry, we now turn to the properties of the radial 
flow considered in this work. As described in the introduction, we model the flow by a relativistic perfect fluid, 
thereby neglecting effects related to viscosity or heat transport, and further assume that the fluid’s energy density 
is sufficiently small so that its self-gravity can be neglected. We also assume that the fluid is minimally coupled to 
gravity, so that it satisfies the standard relativistic Euler equations. 1 The fluid is described by a four-velocity field 
u = u^dfj,, normalized such that = — c 2 , and the particle density ra, energy density e, pressure p , temperature 

T, and entropy per particle s measured by comoving observers. For the following, the enthalpy per particle, defined 
as h := (e + p)/n, will play an important role. We require local thermodynamic equilibrium and the existence of an 
equation of state h = h(p, s) satisfying the first law of thermodynamics 

dh = Tds + (4) 

n 


1 See, for instance, Ref. El for modified gravity theories with non-minimally coupled perfect fluids, in which case a term involving the 
gradient of the Ricci scalar would appear on the right-hand side of Eq. |6|). 
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The dynamics of the fluid flow is governed by the equations 


V M J^ = 0, 
= 0 . 


(5) 

( 6 ) 


wherein is the particle current density and = nhu^u v + pg^ is the stress-energy tensor, and V refers 

to the covariant derivative with respect to the spacetime metric g. Eq. ([5]), the component of Eq. ^ along u and the 
first law, Eq. Q, yield 

Tu^V^s = 0, (7) 

implying that there is no heat transfer between the different fluid elements. 

For the particular case of a radial flow, the four-velocity u has vanishing angular components, so that u = v^dt + 
u r d r = u°k + u r d ri with u° = cu *. The four-velocity’s radial component u r = dr( u) is negative in the accreting case. 
The radial component v of the fluid’s three-velocity measured by static observers is defined by 


1 


u = 


V 1 -1 


(ce 0 +nei), /3 := - 


( 8 ) 


where e 0 := k/ydx is the static observer’s four-velocity and ex := \fNd r is a unit radial vector field orthogonal to it. 
Consequently, the two velocities u r and v are related to each other via 


/ = s/N^j 




(9) 


While u r is well-defined everywhere, the velocity v is only defined outside the horizon, where N(r ) is positive, since 
there are no static observers inside the black hole. A third velocity which will play a very important role in the 
description of the fluid flow is the fluid’s local speed of sound v s , defined as 


dp 

~de 


,9 log h 


9 log n 


( 10 ) 


After these remarks, we derive the equations of motion describing a steady-state radial flow on the spherically 
symmetric static black hole background described by Eq. ([!]). Since the fluid is steady-state and radial, the quantities 
u r , n, h, v s , p, s describing the fluid depend on the areal radius coordinate r only. Eq. 0 immediately implies that 
s = const ; hence each fluid element has the same entropy and it is sufficient to consider an equation of state of the 
form h(n). Next, the particle conservation law, Eq. ([5]), yields 


— (ar 2 nu r ) = 0, 
dr x y 

from which it follows that the particle flux through a sphere of constant areal radius r, 

j n := 47t ar 2 nu r = const., 


( 11 ) 


( 12 ) 


is independent of r. (Notice that j n < 0 is negative for accretion.) Another conservation law is obtained by contracting 
Eq. § with k v and using the Killing equation V M /c„ + = 0. This yields 


d_ 

dr 


-r- - V t = 0 


(13) 


which implies the conservation of the energy flux j e := —4nar 2 T r t /c through the spheres of constant r. Using 
T r t = nhu r ut and the normalization condition of u, we obtain 


j E = 4:Tra 2 r 2 nhu r i 


' N + ( — ) = const. 


(14) 


Summarizing, 2 the equations of motion describing a radial, steady-state perfect fluid flow on a static, spherically 
symmetric background consist of the two conservation laws ( 12|14 ), together with an equation of state h = h(n). 
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Note that it is sufficient to consider Eq. JH| and the components of Eq. ([gJ along u and k, since the angular components of Eq. 
trivial as a consequence of spherical symmetry. 


§ 


are 
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Taking the quotient of Eqs. (12 1 and (14), and eliminating u r from Eq. (12) leads to the following implicit problem: 

,2 


F(r, n ) := h(ny 


cr(r) 




Je 


Jr 


= ( — ) = const., fi:= —, 


(15) 


for the particle density n, where we recall that the function cr(r) = a(r) 2 N(r ) is given by the metric fields and is 
subject to the assumptions (M1)-(M4). Therefore, the problem of determining the flow is reduced to the problem of 
finding an appropriate level curve of the function F which associates to each value of r a unique value of the particle 
density n(r). Once n(r) has been determined, the radial velocity u r is obtained from Eq. (12), 


» 




a(r)r 2 n(r) ’ 


(16) 


from which 


v(r) /x 

c yj y 2 + r 4 <r(r)n(r ) 2 


(17) 


C. Assumptions on the fluid equation of state and main results 


So far, the function h(n), describing the fluid equation of state, has been left unspecified. We now formulate our 
conditions on h and state our main results concerning the solutions of Eq. (15). Their proofs are given in the following 
section. 

We assume h : (0, oo) -A (0, oo) to be a smooth function satisfying the following properties: 


(FI) Positive rest energy: 

As n —> 0, h(n) —> eo converges to a positive number eg > 0. 

(F2) Positive and subluminal sound speed: 

0 < v s (n) < c for all n > 0 . 

(F3) Technical restriction on the derivative of v s [n ): 

0 < W(n) := fgpf < | for all n > 0. 

While the first two conditions are physically well-motivated, the meaning of the third condition is not immediately 
obvious. As we will see in Secs. [TTT] and |IV| it is required to ensure the existence of a unique critical point in the phase 
flow. For the particular case of a polytropic equation of state, 


p = kn 1 , 


(18) 


with k a positive constant and 7 the adiabatic index, integration of the first law, Eq. Q, with ds = 0 gives 


h{n) = e 0 + 1 

7-1 


(19) 


with eo the rest energy of the particle. It is easily verified that the conditions (F1)-(F3) are satisfied in this case, 
provided eo > 0 and 1 < 7 < 5/3. However, more general equations of state with variable adiabatic index 7 are also 
contained in our assumptions (F1)-(F3). 


Before we state our main results concerning the solution of Eq. (15), let us make the following observation: assume 


there exists a smooth solution n(r) of Eq. (15) which is well-defined on the interval [ru, 00 ) including the event horizon 
at r = vh and which converges to a finite, positive value rioo at infinity. By differentiating F(r,n(r)) = const, with 
respect to r we obtain 


n(r)) + ^(r, n{r))n'{r) = 0 . 

The implicit function theorem guarantees the local existence and uniqueness of n(r) as long as 


8F , 1 

ad r>n) = 


2/l(7 


u s / \ / , v s 

- 1 - 72 


( 20 ) 






















6 


is different from zero. According to our assumptions, the right-hand side is negative at the horizon (where a = 0) and 
positive for large r (where a is close to one). Therefore, while the hypotheses of the implicit function theorem are 
satisfied close to the horizon and in the asymptotic region, they are violated at some intermediate point r = r c , where 
dF/dn(r c ,n c ) = 0, n c := n{r c ). In view of Eq. (20) the only possibility for n(r) to have a well-defined first derivative 
at r = r c is that the partial derivative of F with respect to r also vanishes at ( r c ,n c ). In other words, a necessary 
condition for obtaining a differentiable level curve n(r) of Eq. (15) which extends from the horizon to infinity and has 
a finite, positive limit as r —> oo is that this curve passes through a critical point of F. 

Our first main result is concerned with the existence and uniqueness of the critical point (r c ,n c ) of F and the 
structure of the level curves passing through it. 

Theorem 1 (Critical point and level curves passing through it) Under our assumptions (Ml)-(Mf) on the 
background metric and (FI)-(F3) on the fluid equation of state, the function F : [rn, oo) x (0, oo) — > R, (r, n) K > F(r , n) 
defined by Eq. possesses a critical point ( r c ,n c ) for each large enough value of \p\. The critical point (when it 
exists) is unique and describes a saddle point of F at which two level curves of F cross. These level curves can 
be parametrized by (r,n+(r)) and (r,n_(r)), r > th, respectively, with the functions n± satisfying the following 
properties: 

(a) n+ : (r#,oo) —> R is the unique differentiable function satisfying n+(r c ) = n c and F(r,n+(r )) = F(r c ,n c ) for 
all r > th which has regular limits as r —» r B a,nd r —> oo. Moreover, n ^ := linv^oo n+(r) > 0. 

(b) n_ : (rjj, oo) — > R is the unique differentiable function satisfying n_(r c ) = n c and F(r,n_(r )) = F(r c ,n c ) for 
all r > th such that lim,,-^ n_(r) = 0. 

Both functions n± are monotonously decreasing. 


The function n + describes the Michel flow, a radial steady-state accretion flow which extends from the event 
horizon of the black hole to the asymptotic region and which has vanishing radial velocity and finite, positive particle 
density rioo at infinity. An example plot showing the contours of the function F and the level curves (r, n±(r)) for 
the particular case of a polytrope on a Schwarzschilcl background are shown in Fig. [2] below. As will become clear 
from the proof in the next section, the flow is subsonic for r > r c and supersonic for r < r c , the critical radius r c 
representing the sonic sphere. In contrast to this, as we will see, the function n_ diverges as r —» r# and thus it is not 
well-defined at the event horizon. The corresponding flow is subsonic for m < r < r c , supersonic for r > r c , and at 
infinity the particle density vanishes and the radial speed is positive. When the black hole is replaced by a compact, 
spherical star this solution plays an important role in the description of relativistic stellar winds [16] . 

The functions n± in Theorem |T] parametrizing the level set F[r,n ) = F(r c ,n c ) corresponding to the critical point 
(■ r c ,n c ) depend on the accretion rate parameter p, see Eq. (15). Instead of parametrizing the Michel flow by p one 
can also parametrize it in terms of its particle density at infinity rioo- This is the statement of our second main result: 


Theorem 2 Let the conditions (Ml)-(Mf) on the background metric and (F1)-(F3) on the fluid equation of state be 
satisfied. Then, given any n 00 > 0, there exists a unique value of \p\ such that the function F has a critical point with 
the property that the function n+ defined in Theorem [7] satisfies linv^oo n+ (r) = n^. 

The value of n x at infinity can be determined by the particle density of the ambient matter far from the black 
hole. For the particular scenario of a black hole accreting gas from the interstellar medium, it is common to fit the 
flow’s temperature T 00 to the temperature of the interstellar medium, which in turn is related to the sound speed at 
infinity Uoo via the equation of state. For nonrelativistic baryons, for example, we have [lj 

nksT = p = kn 5 / 3 , 

which, taking into account the proton rest mass m p = 938 MeV/c 2 , leads to 


'5k B T a 

Vcv = | — - 


3 


1/2 / T \ 1/2 km 

~ 11.7 f ) -, 

v 10 4 /\ ) s ’ 


which is considerably smaller than the speed of light for typical temperatures of T x — 10 4 Ab 
When Voo/c C 1 we expect that the sonic sphere is located far from the event horizon, r c 
polytropic equation of state with adiabatic index 1 < 7 < 5/3, we confirm this expectation in Sec. IV where we 


rjj. Assuming a 
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compute the leading order contributions in Vao/c to the parameters characterizing the flow. For the particular case 
7 = 5/3 mentioned above we find the following expressions for the compression rates to leading order in Vao/c: 


n c 

rioo 

riH 

rioo 


T (t) 


- 3/2 


~ 2.224 x 10 6 r 3 / 2 (- r ^- S ] 
Vl0 4 A7 


3 / 2 

zho f —)~ 3 ~ 1.334 x IQ 12 ( ZH0 ) 
\cJ V 0.14677/ 


- 3/4 

1 

( 21 ) 

/ rp \ — 3/2 

(io 4 A") 

( 22 ) 


where the dimensionless quantities £ and zho are defined in Eqs. (88) and (90) below. They depend on the parameters 
cri = 2 m/rjj and 02 hr the asymptotic expansion 


r (r) = 1 — 07 —— <72 + O I —rr 


of the metric function a. For a Schwarzschilcl black hole, l = 1 and zho — 0.14677, and the compression rates in 
Eqs. fl2l|22 | reduce to the known results in the literature BE]; however these values may be significantly different 
for distorted black holes. The accretion rate is given to leading order in Voo/c by: 


Je ~ 7T m p c' 




- 3/2 


n 0 


( lCTO 3 / 


J 

5 

s 


(23) 


and only depends on the total mass M of the black hole, not on the details of the metric. 

In the next section, we prove Theorems [l] and [2] by reformulating the implicit problem in Eq. (151 as a (fictitious) 
Hamiltonian dynamical system on the phase space (r, n) . The proof generalizes our previous work m (see also 
Refs. [T3J US)) f° r the particular case of a polytropic fluid accreting on a Schwarzschild black hole background. A 
reader who is not interested in the proof may skip this section and proceed to Sec. m where we discuss specific 
examples involving different equations of states and background black holes. 


III. PHASE SPACE ANALYSIS 


The main objective of this section is to prove Theorems [l] and [2] In a first step, we rewrite the fundamental 
equation (15) in terms of dimensionless quantities and state some preliminary results regarding the behavior of the 
metric fields and the specific enthalpy function in Sec. Next, in Sec. |III B| we reformulate the problem in terms of 

a dynamical system on the phase space (r, n), where the vector field describing the dynamics is the Hamiltonian vector 
field associated with the function F defined in Eq. (15). By construction, F is constant along the phase trajectories 
and thus the latter describe its level curves. The main advantage of casting the problem into a dynamical system 
relies on the fact that the behavior near the critical point of F (where the hypotheses of the implicit function theorem 
break down) can be analyzed using standard tools from the theory of dynamical systems [191 2U:. Using local analysis, 
we show in Sec. IIIB that for large enough values of |/i| the Hamiltonian flow possesses a unique hyperbolic critical 


point with associated stable and unstable local one-dimensional manifolds, which are defined as the set of points lying 
in a vicinity of the critical point which converge to it in the future and past, respectively, along the phase flow. In 


Sec. Ill C we examine the extensions of these manifolds and prove that the unstable branch extends all the way from 
the horizon to the asymptotic region, and thus it describes the physical relevant solution. Our analysis in Sec. |III C| 
also reveals the global structure of the phase space, and in particular proves Theorem [l] Finally, the statement of 
Theorem [2] i s proved in Sec. |III D| The proofs of some technical Lemmata needed for our analysis can be found in 
Appendix |B| 


A. Dimensionless quantities and preliminary results 


We introduce the following dimensionless quantities: 


x := 


r 

th' 



f(z) 


h[n) 

eo 


v\z) 


n _ jn 

. , fl . 2 

no ATTcrjjTiQ 

9 log f{z) <9 log 7z) 

— X~j -, w{z) := —=:- 

a log z o log z 


(24) 

(25) 




























Here, no denotes a typical density (which could be defined as the particle density at infinity, for example), fi is the 
dimensionless accretion rate, v is the local speed of sound divided by c, and / denotes the specific enthalpy. Note 
that in the accreting case, both fi and u are negative. For notational simplicity we shall drop the bar on fi in the 
following. In terms of these quantities, the fundamental equation, Eq. (151, can be rewritten as 


Fn(x,z) = 


Je 

jn^ 0 


= const., 


(26) 


where F ^ : f 2 


is the function on fl := ( 0 , oo) x ( 0 , oo) defined by 


F^(x,z) := f(z) 2 


cr(x ) + 




(x, z) € fl. 


(27) 


For the analysis of the Hamiltonian flow associated with the function F^ we will need two technical lemmata. The 
first one is related to the behavior of the metric field a. 


Lemma 1 The function 


i(x) := 


r(x) 


xcr'(x) 

is a monotonously increasing function satisfying a(l) = 0 , lim. 


0 < a(x) < r-^-x, 
2m 


x —^OO 
2 


x > 1, 

a(x)/x = rn/i^rn) and the inequalities 


(4a + 1) 
16a+ 1 


< xo! < 1 + 4a, 


(28) 


for all x > 1. 

Proof. See Appendix [B| p 

The second lemma, which is a consequence of the assumptions (F1)-(F3) on the fluid equation of state, gives 
information about the asymptotic properties of the specific enthalpy / and the dimensionless sound speed v as z —> 0 
and 2 — i oo. 

Lemma 2 (a) Given z\ > 0 there exists S > 0 such that the following inequalities hold for all 0 < z < Z\: 


v{z) > fc 1 / 3 , 


f{z) > exp 


3 

2 


5 2 z 2 / 3 


(29) 


(b) v is a strictly monotonously increasing function satisfying 

lim v(z) = 0 , lim v(z) = V\ 

z—t 0 z—> oo 


for some constant 0 < < 1. 

(c) Given zq > 0 there exists a constant q > 0 such that for all z > zq, 

/w>/m(|) - 

In particular, lim f(z) = oo. 

z—t OO 

Proof. See Appendix [B| 


(30) 


□ 
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B. Definition of the dynamical system and critical points 


The Hamiltonian vector field with respect to the function F / is defined by 


X F (x,z) := (VF Al (x,z))- L 


dF^(x,z) 

dz 

dF^(x,z) 

dx 


(x, Z ) € fi. 


(31) 


The associated Hamiltonian flow is determined by the integral curves of Xp , obtained by solving the ordinary differ¬ 
ential system 


d 

dX 



X F (x,z) = f{z) 2 




(32) 


with A a fictitious time parameter. By construction, these curves leave the level sets of F t , invariant, so that the 
problem of finding implicit solutions of the equation F ;/ (x, z) = const, can be understood by determining the qualitative 
properties of the phase flow. 

The critical points (x c , z c ) of the system are determined by the points at which X F vanishes, that is 


v(z) 2 (cr(x) + u 2 ) — u 2 = 0, 
4 y 2 

a'{x) - — = 0, 
x 


(33) 

(34) 


where u 2 = /r 2 /(x 4 2 2 ). Eliminating z, solving this system is equivalent to finding the zeros of the auxiliary function 


J" M (x) := v(z(x)) 2 [1 + 4a(x)] - 1, 


z{x) 


M 

\Jx b a'{x) ’ 


x > 1, 


(35) 


where a(x) = a(x)/(xa 1 (x)) is the function defined in Lemma [T] We first show the existence of a zero by analyzing 
the limits x —> 1 and x —> oo. For x —> 1, 


^(1) = ^(2| A *|/ X /^(I)) 2 -1<0, 


where we have used the assumptions (M2) and (F2). In order to analyze the asymptotic behavior of X^x) for 
large x, we first note that for any fixed value of fj,, z(x) —> 0 as x —> oo since according to assumption (Ml), 
x 2 cj'{x) —> cr-| = 2m/r F ■ Next, we fix the constants z\ and 6 in Lemmata) which only depend on the equation of 
state. For any fixed value of n and large enough x we then have 


Xfj,(x ) > S 2 z(x) 2 ^ 3 [1 + 4a(x)] — 1 = S 2 


(2|^|) 2/3 

(a^er^x )) 1 / 3 


1 

X 


4cr(x) 

2 <t'(x) 



- 1 . 


Therefore, we can guarantee that X^(x) is positive for large x if we choose |/i| large enough such that |/i| > af/(16S 3 ). 
In this case, the smooth function T tl is negative for x = 1 and positive for large x, and thus it must have a zero by 
the intermediate value theorem. 

Next, we show that this zero must be unique. For this, we consider the derivative of 


X/(x) = v(z(x)) 2 


w(z(x)) 

xa(x) 


[1 + 4a(x)][xa'(x) — 1 — 4a(x)] + 4a'(x) 


(36) 


It follows from our assumptions (M4) and (F3) that the right-hand side is always positive. Indeed, when w = 0 the 
expression inside the curly parenthesis reduces to 4a'(x) > 0, while for the other extreme, w = 1/3, it reduces to 

-[(1 + 4 a){xa! — 1 — 4a) + 12xaa , l = - T(1 + 16a)xa / — (1 + 4a) 2 l , 

3xa 3xa 

which is also positive as a consequence of Lemma |TJ Therefore, the function is strictly monotonously increasing, 
and as a consequence, the critical points of the dynamical system are unique. 

After having determined the critical point we analyze the flow in the vicinity of it. To this purpose, we linearize 
the Hamiltonian vector field X F at the critical point (x c , z c ). The linearization is given by the matrix 


DX f (x c , Z(dj 


d 2 Fn(x 0 ,z 0 ) d 2 F,(i c ,z c ) 

d 2 F fi (x c ,z c ) d 2 F^(x c ,z c ) 
dx 2 dxdz 


cr'(x c ) 


f(Zc ) 2 

Z c 


2 

[x c a'{x c ) - 1 - 4a(x c )] 



(37) 
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where A := zjx c and r/ := 1 — u(z c ) 2 +w(z c ) > 0. A short calculation which exploits the fact that v(z c ) 2 [l+4a(x c )] = 1 
reveals that the determinant is 


det [DX f (x c ,z c )\ = 


-x c a'{x c ) 


2 /Oc) 4 




)< 0 . 


(38) 


Since DXp(x c , z c ) has zero trace, it follows that its eigenvalues are real and have different signs; in particular ( x c , z c ) 
is hyperbolic and the Hartman-Grobman linearization theorem m is applicable. Furthermore, the negative and 
positive eigenvalues of DxF(x c , z c ) give rise to one-dimensional stable and unstable local manifolds, corresponding, 
respectively, to the set of points (x,z) converging to the critical point when A —> oo and A —> — oo. Accordingly, the 
function A), has a saddle point at ( x c ,z c ), since the determinant of DXp(x c , z c ) is equal to the one of the Hessian 
of F,j at ( x c ,z c ). In a vicinity of (x c ,z c ), the critical point together with the stable and unstable manifolds describe 
the level set of F ^ through the critical point. In the following, we study the extensions of the stable and unstable 
manifolds to the whole phase space H. 


C. Extension of the stable and unstable local manifolds and global properties of the phase space 


In the previous subsection, we showed that for large enough |/i| there exists a critical point of the dynamical 
system (32), and further we showed that a critical point is always unique and hyperbolic, with corresponding one- 
dimensional stable and unstable local manifolds T_ and F + . In this subsection, we discuss the extensions of these 
manifolds towards the horizon (x = 1) and towards the asymptotic region x —> oo, and we discuss the global qualitative 
behavior of the phase flow. 

Our result is based on the analysis of the two curves Ti and r 2 , corresponding to the points ( x,z ) with vanishing 
partial derivative of F^(x,z) with respect to x and 2 , respectively. Using Eqs. ( 33|34 ), these two curves are therefore 
parametrized by 


zi(x) := 


2H 


\J x^ct' (x) 

with the function H : (0, oo) —> (0, oo) defined as 


z 2 (x) := H 1 (xct(x)), x > 1, 


H{z) := 


A 


1 


- 1 


Notice that H is invertible since H{z ) 


dH{z) 

dz 


v{z ) 2 

oo as z —> 0 and H(z ) 
2 


, z > 0. 

0 as z —> oo, and since 


(39) 


(40) 


FhFlz) 


[l - V 2 {z) + w(z)] < 0, 


which shows that H(z) is monotonously decreasing. Furthermore, using Lemma [I] 

dz\ . Zi(x) xa'(x) — Aa(x) — 1 

* (x) = -iw- <0 ’ 


dz 2 . _ X 3 


xa' (x) + 4tr(x)] < 0. 


Therefore, ^i(x) and z 2 (x) are monotonously decreasing functions. The curves Ti and T 2 intersect themselves only 
at the unique critical point (x c ,z c ), and therefore, they divide the phase space into four regions, see Fig. [ljbelow. In 
each of these four regions, the flow has a particular direction which is indicated by the arrows in this figure. We note 
in passing that the curve T 2 has a nice physical interpretation: it divides the phase space into regions where the flow 
is sub- and supersonic, respectively. In order to see this, recall the expression ( |17[ ) for the radial velocity of the flow 
measured by a static observer. In terms of dimensionless quantities, 


Comparing this expression with 


v(z 2 (x)) = 


yj H 2 + X A <j'{x)z 2 

_ \jA _ 

yV 2 + x 4 cr'(x)z 2 (x) 2 
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which follows from the definition of Z 2 {x ), and using the monotonicity of v(z), this implies that \v\/c < v for points 
(x,z) lying above the curve 1?2, while \v\/c > v for points (x,z) lying below this curve. Since z^{x) —>■ oo as x —> 1, 
an important conclusions is that any differentiable solution z(x) of the implicit equation (26) which extends from the 
event horizon to x — >• oo must be transonic. As discussed in Sec. [TT] this solution must cross the curve T 2 precisely at 
the critical point, and hence it must coincide with either the stable or unstable local manifold. 



FIG. 1: A sketch of the phase space, showing the critical point ( x c , z c ), the stable (r_) and unstable (r+) manifolds, and the 
curves Fi and F 2 which delimit the four regions I-IV. Also shown is the direction of the flow in each of the regions and along 
the curves Ti and IV The flow’s radial velocity measured by static observers is subsonic in the region above the curve F 2 and 
supersonic in the region below F 2 . 


In order to proceed we compare the locations of the stable (T_) and unstable (T + ) local manifolds with those of 
the curves Fi and in the phase space Q. For this, we first compute the slopes of these four curves at the critical 
point. The slopes of r± at the critical point are determined by the eigenvectors of the linearization DXp(x c , z c ). 
From Eq. (37) one finds the eigenvectors 


e± = 


-^(1 Tk ) 


corresponding to the eigenvalues ±2a'(x c )f(z c ) 2 k/z c , where 


k= \j 1 + [x c a'{x c ) - 1 — 4a(x c )]. 

Since DXp(x c , z c ) has negative determinant, k is real and we can take the positive root, and according to Lemmajl] 
k < 1. Therefore, the slopes of the stable and unstable local manifolds are given by 


—(1 + fc), 4 = - — (1 - k), 
r] V 

respectively. The slopes of Ti and T 2 can be obtained by direct computation; the result is 

h := d ^(x c ) = - — (l-k)(l + k), 
ax rj 


(41) 


(42) 


; dz2 r \ 
h ■■= ~r{Xc) = ^ 
ax 


and 


2A 

V 


(43) 
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Eqs. ( 41 | 42 | 43 ) and the fact that 0 < k < 1 imply the inequalities 


l- <l 2 <h < 1 + < 0 . ( 44 ) 

Therefore, close to the critical point, the relative location of Ti, T2 , r+ and T_ is as shown in Fig. [Tj In particular, 
the curves Ti and T 2 lie between the stable and unstable local manifolds, at least in a vicinity of the critical point. 
In the following, we show that this is also true globally, that is, T± do not intersect the curves Id and T 2 at points 
other than the critical one. In order to prove this we make the following general observations: 


(I) Consider region I, which is the region that lies above the curves Ti and T 2 . Here, > 0 and ^7 > 0 , implying 
that x increases and 2 decreases along the flow of the dynamical system ( 32 ). Moreover, at the boundary points 
x > x c , 2 = zi(x), the flow is pointing inward , implying that it cannot leave region I at those points. As a 
consequence, the unstable manifold T + must extend to x —>■ 00, since according to Eq. ( 31 ) x increases and z 
decreases but is bounded from below by 2 > z±(x). 


Using similar arguments for the time-reversed flow, one concludes that the stable manifold T_ extends from 
(x c , z c ) to 2 —> 00 inside region I, with x decreasing monotonically. 


(II) Next, consider region II which lies below the curves Ti and T 2 . Here, ^ < 0 and §7 < 0 , implying that x 
decreases and 2 increases along the flow of the dynamical system ( 32 ). At the boundary points x < x c , z = z\{x) 
the flow is pointing inward , implying that T + cannot leave region II at those points. Since x decreases and 2 
increases, and since the curve Ti crosses the horizon at x = 1, it follows that r + extends from (x e , z c ) to a point 
(1, Zh) where it crosses the horizon, with z c < Zh < 21 (1). 


By a similar argument, one can show that the stable manifold T_ extends in the past from (x c , z c ) to x —> 00, 
2 —> 0, since it must lie between the x-axis and the curve T 2 . 


In conclusion, we have shown that the unstable manifold extends from the the asymptotic region x —> 00 through 
the critical point to the horizon, where 2 ( 1 ) = Zh has a finite value. The corresponding flow is subsonic for x > x c 
and supersonic for x < x c . In contrast, the stable manifold extends from x —> 00 through the critical point, but z(x) 
diverges as x — > 1 approaches the horizon. The corresponding flow is supersonic for x > x c and subsonic for x < x c . 
Other possible flows are given by integral curves which lie below the curves T + and T_, whose corresponding flow is 
everywhere supersonic, and integral curves which lie above the curves T + and T_, whose flow is everywhere subsonic. 
In all cases, 2 increases as x decreases and hence, one can associate to each value of the dimensionless radius x = r/rn 
a unique value n(r) = noz(x ) of the particle density, which is a monotonously decreasing function of r. 

In order to conclude the proof of Theorem [l] it remains to analyze the asymptotic behavior of the functions z ±, 
parametrizing the manifolds T±, for x — > 00. Let us start with 2+. Since 2+ is monotonously decreasing and since 
2-|_ (x) > 2i(x) for all x > x c the limit 


Zoo 


lim 2_|_(x) 


exists. As we argue now, this limit cannot be zero. Indeed, 2+(x) > z\ (x) for all x > x c implies that along r + the 
radial velocity is bounded by 

I, _M_ M = Vx 2 a'(x) 

1 1 X 2 2 + x 2 2i(x) 2 y/x 

in the asymptotic region. According to assumption (Ml), the numerator converges to a positive constant for x 00, 
implying that u — > 0 as x —> 00. Therefore, since is constant along T + , 

Fn(xc, z c ) = F^Xc, 2+(x c )) = lim F^(x, 2+(x)) = /(2oo) 2 , (45) 

x—too 

where we have used Eq. ( | 27 [ ) in the last step. Lemma below shows that i^(x c ,2 c ) > 1 , implying that 200 must 
be positive. Therefore, along T + the particle density n[r) = noz(x) converges to a finite, positive value and the 
radial component u r of the four-velocity and the fluid’s velocity v measured by static observers converge to zero when 
r —i 00 , as expected. 

In contrast to this, along the stable manifold T_ we must have 


lim 2_(x) = 0, 

X—>00 


since 0 < 2_(x) < 2 2 (x) —i 0 as x —► 00, implying that 

1 < F m (x c , 2 c ) = lim F /1 (x,2_(x)) = l+ulc, 

x—too 

and while the particle density n{r) = n^z^x) converges to zero, u r and v converge to non-zero values as r — > 00. This 
concludes the proof of Theorem [l] 
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D. The matching problem 


In the previous subsections, we discussed the local and global properties of the flow for a given value of the 
dimensionless accretion rate p. In this subsection, we prove Theorem [2] which states that for any given positive value 
rioo > 0 of the particle density at infinity there exists a unique value of p such that the function F ^ has a critical 
point whose function n+ parametrizing the unstable manifold satisfies 

lim 7i+(r) = noo. 


So let Zoo > 0. According to Eq. (45) we need to find p and a critical point ( x c , z c ) of F ^ such that 

= f(Zoof, 


E^(x c , z c ) — f (zf) 


2 Cr (2 ; c) _ \ 2 


i - v(z c y 


where we have used Eq. (J33J) in the first equ ality. Instead of p, we parametrize the critical point by z c and use the 
fact that v(z c ) 2 [ 1 + 4 a(a; c )]== 1 (see Eq. (35)) in order to express z c as a function of x c . Therefore, we introduce the 
auxiliary function C : (0, oo) —R defined as 


C(z) ~ f{z) 


2 <r{x c (z)) 
l-viz) 2 ' 


z > 0, 


with the function x c : (0, oo) —► (1, oo) given by 

x c (z) := a -1 Q H^T 2 - l]^ , z>0. 


(46) 


(47) 


where the inverse of the function a : (1, oo) —> (0, oo) exists according to Lemma]!] Our goal is to show the existence 
of a unique z c > 0 such that 


C{Z C ) = F^XcZc) = /(Zoo) 2 - 


(48) 


According to the result of Lemmata), the right-hand side is strictly larger than one. The result then follows from: 
Lemma 3 The function C : (0, oo) —» R defined in Eq. (46) is strictly monotonously increasing and satisfies 

lim C(z) = oo 


lim£(z) = 1, 

z—>• 0 


(49) 


f 2 ■ 

> OO 


In particular, there exists for any /oo > 1 a unique z c > 0 such that C{z c ) = / c 

Proof. See Appendix [B] p 

As a consequence of Lemma [3j given z^ > 0, there exists a unique z c > 0 such that C(z c ) = /(zoo) 2 . Defi ning 
(x c ,z c ) := (, x c (z c ),z c ) and p := ±^x^.a'(x c )z c / 2 it is not difficult to verify that the function F ^ defined in Eq. (27) 
has a critical point at (x c ,z c ). Since |/i| depends uniquely on x c and z c , and x c depends uniquely on z c , ji\ is 
unique. According to Theorem [l] there exists a unique differentiable function z+ : (1, oo) —► R such that z+(x c ) = z c , 
F^(x, z+(x)) = F^Xc, z c ) for all x > 1 and which has regular limits as x —> 1 and x —> oo. Since 

/ ( lim z+(a;)) = lim F^x, z + (x)) = Fjx c ,z c ) = £(z c ) = f(z oo) 2 , 

\x—too / x —>-oo 

it follows by the monotonicity of / that 


lim z+(x) = z c 

x—yoo 


as desired. 
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IV. EXAMPLES AND COUNTER-EXAMPLES 


In this section, we discuss several particular examples of accretion flows. We start with the study of the two limiting 
cases in condition (F2), corresponding to perfect fluids with sound speed v s = 0 and v s = c, respectively. In the first 
case, h = e o is constant implying that the gas does not have internal energy nor pressure, and the matter in this 
case is dust. The second case describes a stiff fluid, for which the enthalpy is proportional to the particle density, 
h(n) = kn for some positive constant k. Although in both cases the function F defined in Eq. (15) does not possess 
any critical point, there are still non-trivial, globally defined solutions which we discuss below. Next, we consider 
three numerical examples concerning the radial accretion of a polytropic fluid. The first two examples are for a 
Schwarzschild background and have adiabatic index 7 = 4/3 and 7 = 1.99, respectively. In the first case, one obtains 
the standard Michel solution while in the second case which violates the condition (F3) there is no global solution. 
The third example has adiabatic index 7 = 5/3 and describes a flow on a distorted black hole which violates condition 
(M4). While there is still a global flow, it has features which are different than in the standard case where all the 
assumptions are satisfied. Finally, we analyze a radial polytropic fluid flow with 1 < 7 < 5/3 on an arbitrary static, 
spherically symmetric black hole satisfying the conditions (M1)-(M4) assuming that the sound speed at infinity i>oo 
is much less than the speed of light. In this case, the sonic sphere is located far from the event horizon, r c r#, and 
the parameters characterizing the flow, like its accretion rate, compression rates etc. can be expanded in powers of 
Voo/C. 


A. Accretion of dust 


For dust, h(n) = eo is equal to the rest energy of the particles, and in this case Eqs. (151 and yield the simple 
relation 


cr(r) 


1 - 


v(r) 


Js 


= const. 


(50) 


Taking the limit r —> 00 and requiring that v{r) —> 0 and <r(r) —> 1 in this limit, the constant must be equal to one 
and it follows that the energy and particle fluxes are related to each other by j E = eoj n . Then, Eq. (50) yields the 
following expression for the fluid’s velocity measured by static observers: 


j(r) = -csj 1 - cr(r), 


(51) 


or v{r) = —\/2GM/r for a Schwarzschild black hole of mass M, which is just the Newtonian expression for a freely 
falling radial particle with zero velocity at infinity. The particle density can be obtained using Eq. ©> 


n(r) = 


\jn 1 1 

47tc r 2 ^/1 — <r(r) 


(52) 


or n(r) = \j n |/(47 tV 2 GMr 3 ) in the Schwarzschild case. The condition (Ml) for the existence of a positive mass at 
infinity implies that n(r) — > 0 as r —> 00, so the particle density at infinity is zero. In fact, there are no solutions 
of Eq. (50) with > 0: the assumption rioo > 0 leads to je/i^ojn) = 1 when taking the limit r —> 00 in Eq. (15), 


and it follows that n(r) is again of the form Eq. (52), leading to a contradiction. Therefore, the existence of a radial 


steady-state fluid accretion flow into a nonrotating black hole requires nonzero pressure for the particle density at 
infinity to be positive. 


We also see from Eq. (52) that there is no global regular solution if a is not restricted to be everywhere strictly 
smaller than one. 


B. Accretion of a stiff fluid 


In the stiff case, h(n) = kn for some positive constant k and Eq. (15) yields 

2 / • \ 2 

n 2 cr{r) + —. = ( ) = const , 

r 4 \kjnj 


(53) 
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where /r = j n /(4:Trc). Taking the limit r —» oo and assuming again that er(r) —> 1 in this limit gives rioo = j e /{kj n ), 
which hxes the relation between the energy and particle fluxes. Next, taking the limit r —> th in Eq. (53) yields 
H = —r^rioo, which determines the particle and energy fluxes to be 


j n = -Aircrnrioo, j e = -^ck^Hn^) 2 . 


With these constants Eq. (53) yields the expression 


i{r) = 


' 1 - 


r(r) 


for the particle density. Using l’Hopital’s rule, we obtain the compression rate 

n(r H ) 2 


rioo \/ r H a'{r H ) 


(54) 


(55) 


(56) 


between the particle densities at the horizon and at infinity, which is finite if the condition (M2) is imposed. For the 
Schwarzschild case, 


i(r) = ■ 




4m 2 ' 




(57) 


and the compression rate is 2. However, it could be considerably lower or higher for black holes exhibiting a different 
behavior close to the event horizon. 


Finally, the velocity measured by static observers is obtained from Eqs. (17) and (551: 

v{r) = -c-f, 


(58) 


which decays as 7 ' 2 (as opposed to r x ! 2 in the dust case). 


C. Polytropic fluids 


Next, we consider a polytropic equation of state, for which 

h(n) = e 0 + :n 7_1 (59) 

7-1 

with k > 0 a positive constant, eg the rest energy of the particle and 7 > 1 the adiabatic index. Condition (FI) is 
satisfied as long as eg > 0. The sound speed is 


s(n) = cJ d l ° gh = c^/ 7 — 1 ,/- 
I/ dlog n \ 7 


nn- 1 


(-V— 1 5 


(60) 


which satisfies condition (F2) as long as 1 < 7 < 2. The logarithmic derivative of the sound speed is 

d log v s 7 — 1 1 


W(n) = 


9 log n 


2 1 + 


7 k 


-rft- 


(7—l) e o 

which satisfies the required bound in (F3), provided that 1 < 7 < 5/3. Defining 

'(7-l)eo ll/(T_11 


n 0 := 


7 k 


we obtain, in terms of dimensionless quantities, z = n/ng, and 


v{z) 

w(z) 


h(n ) 
eg 


1 + z 


7-1 


?£M = v 7vt 

c 

7-1 

2(1 + z'y~ 1 )' 


Z 7-i 

1 + z'y - 1 ’ 


(61) 


(62) 


(63) 

(64) 


(65) 
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We provide three examples showing the level sets of the function F / defined in Eq. (27). The first two examples are 
shown in Fig. [2] and refer to a Scliwarzschild black hole background. In the first example, /i = —2 and 7 = 4/3, and 
there is a unique hyperbolic critical point which is located at (x c , z c ) = (2.20,1.23) through which two level curves 
of F fl cross each other. Only one of these two curves extends from the horizon to infinity, and describes the physical 
accretion flow. In the second example, /r = —0.33 and 7 = 1.99, which violates the upper bound in condition (F3). 
Here, there are two critical points, one corresponding to a saddle point and the other to a local extremum of F tl . 
The right panel of Fig. [2] suggests that the level curves of I7 through the saddle point connect to each other instead 
of extending to the asymptotic region x —> 00. Consequently, there is no global flow in this case. Recently, such 
“liomoclinic-type” accretion flow solutions have also been found in a cosmological context, see Refs. Ennui- 
The third example, shown in Fig. [3j has parameter values /i = —1 and 7 = 5/3 and has a background metric 
described by a distorted black hole, such that 


. . 1 1 4 8 4 

a(x) = \- - - — + — - — 


x > 1. 


( 66 ) 


In this case, the conditions (Ml)-(M3) on the metric are satisfied, however (M4) is violated. An analysis of the 
function FAx) defined in Eq. (35) reveals the existence of three critical points, two corresponding to a saddle point 


of F tJ and one to a local extremum. Unlike the previous example with 7 = 2 there does exists a global accretion flow; 
however, due to the presence of the extremum, the features of this flow are different from the standard situation where 
all assumptions (M1)-(M4) are satisfies. For example, the particle density ceases to be a monotonously decreasing 
function of r. 




FIG. 2: The level sets of the function J 7 / for a polytropic fluid accreted by a Schwarzschild black hole. The level set 
corresponding to the critical point is displayed with the thick black line. Left panel: In this case, 7 = 4/3 and there is a unique 
critical point located at ( x c ,z c ) = (2.20,1.23) corresponding to a parameter value oi fi = —2. The thick black curve is the 
union of the level curves (r, n+(r)) and (r, n_(r)) described in Tlieorem[l] and the former curve which extends from the horizon 
(x = 1) to the asymptotic region (x —» + 00 ) describes the Michel flow. Right panel: In this case, 7 = 1.99 and there are two 
critical points for fj, = —0.33, one corresponds to a saddle point of Fp and is located at (x c ,z c ) = (2.50,0.17), while the other 
one is a local extremum of F/ which is located at (x,z) = (3.7,0.09). In this case, the two level curves through the saddle 
point do not extend to infinity. Instead, they seem to wind around the local extremum and connect to each other. There is no 
globally-defined radial flow solution in this case. 


D. Small sound speed expansion 

In this subsection, we consider the situation where the sonic sphere lies in the region where the gravitational field is 
weak, that is, r c tr- Furthermore, we assume that in the asymptotic region, the sound speed Uqo is much smaller 
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FIG. 3: The level sets of the function F\ for a polytropic fluid with adiabatic index 7 = 5/3 and /r = —1 accreted by a strongly 
distorted black hole described by Eq. ( 66 l. In this case, there are three critical points, two of which are located at (x, z) = 
(1.16,2.88) and (x,z) = (3.42,0.24), respectively, and correspond to a saddle point, and one located at ( x,z ) = (1.66,1.19) 
which corresponds to a local extremum of F lt . The level set corresponding to the saddle point closest to the horizon is displayed 
with the thick black line. The branch that extends from the horizon to infinity describes a radial, transonic accretion flow. 
However, as can be seen from the graph, unlike the standard case where the condition (M4) holds, in this case the particle 
density is not monotonously decaying in r, but reaches a maximum value just before the flow passes through the critical point 
and becomes supersonic. 


than the speed of light, such that := Uoo/c -C 1. Finally, we also assume that the fluid is described by a polytropic 
equation of state, as in the previous subsection. With these assumptions, we can expand the metric function 




1 


01 = 


2m 

r H ’ 


(67) 


for x > x c , and also we can expand the fluid quantities in terms of the small parameter v^. Our first goal is to 
express the quantities x c , z c , v c , and ijl characterizing the flow in terms of 

We begin with the relation between the particle density at the critical point z c and the sound speed at infinity. 
This relation is determined by the equations (46|48 1, that is, 


CM = fM ^^ll = /(**,)”, 

1 - v{z c y 


( 68 ) 


with the functions / and v given by Eqs. (63) and (64), and the function x c defined in Eq. (47). Using the expansion 
in Eq. (67) we find the following expansion for x c (z c ): 


01 


2(^0) 


= 4(7 - 1 )zj 1 U - 


37-2 + 8(7 — 1 )^| 


.r i +o( 2c 2 ^- 1 )) j. 


Substituting this result in Eq. (68) we find 


C(z c ) = 1 + (5 - 37 )z2 1 + 


1 + 3(7-1)(37-4)+ 16(7-1) 


2 


'll 


^ 1 ) + o(^- 1 )). 


(69) 


(70) 


On the other hand, from Eq. (64) we obtain 


Zoo — 


7-1 


(7 -1 y 


+ 0 ( 1 4) 


(71) 
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for Voo <C 1 , and hence 


f ( Zoo y = (i + z ^f = [i + 


i - 1 - k 


= 1 + 2 


+ 3 


7-1 (7 ~ 1) : 


+ 0{vl). 


(72) 


Comparing Eq. ( |70| ) with Eq. (72) we conclude that for 1 < 7 < 5/3 the solution z c of the equation £(z c ) = /(+x>) 2 , 
which is unique by the nronotonicity of £, must satisfy z c <l. Depending on whether 1 < 7 < 5 / 3 or 7 = 5/3we 
obtain the following behavior for the quantities z c , x c , v c , zh, z^ and p. in terms of v^ -C 1 : 


Case I (1 < 7 < 5/3): 


= 


x r = 


2vi 


.(5 — 37 ) (7 — 1)_ 
5 — 37 cr 1 


[l + Avlc + O(^)] , 


8 vl 


[1 + B^ + 0(0] 


ImI = 


z H = 


4 A 


v oo [l + C^lo + ®{ U to)\ > 

[l + Du^ + C?(i/£,)] 

1 


5-37 

^5-37 \ V+ 


7-1 

«4 ~ 37 

7-1 


1+lD- 


7-1 


*4+0(4,) 


(73) 

(74) 

(75) 

(76) 

(77) 


with 


9-7-1 B — 3-f 

A = 2 5 Pt=ij (5 - 37 ) 


and where the coefficients A, B, C and D are defined by 
1 


A 

B 

C 

D 


(5-3 7 ) 2 (7-l) 2 

2 


12-i(37+l) 2 -32( 7 -l) 2CT § 
£ (J 1 


(5-3y)(7-l) 


37-2 + 8(7-l)^| 


^ 7 1)A 2 (5 — 3 7 ) 2 


9(7 — 3 7 ) + 32(3 — 7 )—| 


1 1 11 

2 7 j (5-3 7 )( 7 -l)~ 4(5-37 ) 2 


3(37+1) + 64(7 — 1)^| 


a + \b- 


8 0-2 _ 1 1 

5 — 37 of 4 (5 — 37 X 7 — l ) 2 


27 7 2 - 667 + 47 + 64(7 - 1) 


2 a 2 


’1 J 


(78) 

(79) 

(80) 
(81) 


Here, x c := x c (z c ) and v c := v{z c ) have been computed by substituting the expansion for z c into Eq. (69) and 
Eq. (641, respectively. The dimensionless accretion rate is obtained by substituting the expansions for x c and z c 
into the expression 2 \fi\ = y/x%a'(x c )z c . The density at the horizon zh is computed from the relation between 
zh and z^ which follows from 


(1 + zJj 1 ) 2 ^- = F^l, z H ) = f(zoo) 2 = (1 + 4 X ) 2 > 


leading to the equation 


M (1 

ZH 


7— 1 
Z H 


) — l + ^So 1 — 1 + 


7-1 (7-1)' 


+ 0 ( 4 )- 


(82) 


Since the dimensionless accretion rate \^i\ converges to zero as fast as i4o 37 ^^ 7 for v a0 -+ 0, Zh must converge 
to zero at the same rate. With this observation in mind, the expansion of zh in terms of follows easily from 
Eq. |82|. 
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• Case II (7 = 5/3): 


(^) 3/2 [1 + 0(01, 

7 ,° 1 - [1 + O(^oo)] , 

8 Voo 


z H = 



where 


l \l 1 + ( x ) — 2 > 


and zho is uniquely determined by the equation 


ZHO - mZHO + 1 ) = 0 , 


?■■= f 8 


3/2 


whose solution can be written explicitly as 

^2,3 


ZH 0 = [ P + W~\ - 

W 


= /? 1/3 / 3 2 + 7 + V3\ P + 2 


1/3 


(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 


A plot of zho as a function of 07 is shown in Fig. [4j In contrast to case I, here the accretion rate converges to 
a finite, positive value as —> 0 and consequently, it follows from Eq. (821 that zr must also converge to a 

finite positive value zho when —> 0. This value is determined by Eq. (89) which follows from Eq. (82) by 

taking the limit v^ —» 0. In the Schwarzschild case, l = 1 and zr 0 — 0.14677. 



2 m 


FIG. 4: The parameter zho vs. <ti = 2m/rH■ 
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We see that in both cases the critical point is located far from the event horizon, r c th when <C 1, which is 
consistent with our assumption of the critical point lying in the weak field region. Finally, we compute the accretion 
rate j e and the compression rates ndriao and njj/n <*, the gas undergoes when it falls from infinity through the critical 
point into the black hole. These quantities can be obtained by combining the expansions above for (i, z c and zh with 
Eq. |7lj ). The result is 


• Case I (1 < 7 < 5/3): 


je 


z c 

Zoo 

zh 

Zqo 


AXirceg 
2 


GM' 

1 

7—1 


5-3 7 




n 00^00 


1+ [A- 


1 +[D + 


(7-1 ) 5 


7-2 
(7-l) 2 

1& + OM 




1 + ID + 


7-2 
(7- l ) 5 


V 2 oo + O(vi) 


• Case II (7 = 5/3): 


/ KjrlVl \ _q 

j e = 7 rce 0 I —I rioo^oo [I + 0[Voo)\, 

/o\ 3/2 

— = [JJ Zho^II + O^cc)}, 


(91) 

(92) 

(93) 

(94) 

(95) 

(96) 


where we have used 77*07 = 2 GM/c 2 with M the mass of the black hole. We see that in both cases the accretion 
rate scales like M 2 n 00 v and that zh is larger than z x by a factor which is proportional to implying that 
significant compression rates are achieved during the accretion process. On the other hand, in case I the compression 
rate at the sonic sphere is almost independent of for small while it scales like Voo ' in case II. In both cases 
the compression rate at the horizon becomes larger as ay = 2 m/ru increases. 


V. CONCLUSIONS 


In this work, we have studied the Michel flow, describing the steady spherically symmetric accretion of a perfect 
fluid into a nonrotating black hole, under the assumption that the fluid’s density is sufficiently small so that its 
self-gravity can be neglected. Our treatment does not assume a polytropic equation of state nor a Schwarzschild 
background geometry. Instead, we imposed rather general conditions on the fluid equation of state and considered a 
large class of static, spherically symmetric and asymptotically flat black holes with a regular horizon. Consequently, 
our results are also applicable for certain deformed Schwarzschild black holes which might arise in alternative theories 
of gravity or in classical general relativity in the presence of external matter fields such as dark matter, for example. 

Our main result is that for a given equation of state and metric satisfying our assumptions, there exists for each 
positive rioo > 0 a unique, spherical and stationary solution of the general relativistic Euler equations which is 
everywhere regular at and outside the event horizon and whose particle density at infinity is rioo. This solution 
describes a transonic flow, the flow being supersonic close to the black hole and subsonic in the asymptotic region. 
We have also computed the accretion rate and the compression rates the fluid undergoes at the event horizon and the 
sonic sphere, assuming a polytropic equation of state and the usual assumption that the fluid’s speed of sound Voo is 
much smaller than the speed of light at infinity. While to leading order in Vao/c the accretion rate is given by the 
familiar Bondi formula and only depends on the background metric through its total mass, the compression rate at 
the horizon depends on the details of the metric fields describing the black hole; in particular it depends on the ratio 
fs/fH between its Schwarzschild radius r s = 2m and the radius of the event horizon ry. For a Schwarzschild black 
hole, r s = ru and this ratio is one. However, a spherically symmetric static black hole in the presence of external 
fields satisfying the weak energy condition has r s > r # (see Eq. (All), leading to higher compression rates at the 


horizon than in the absence of external fields. Intuitively, this might 
surface gravity, and in this sense the gravitational acceleration at the horizon is larger when r# < 
Schwarzschild case, resulting in higher compression rates. 


re expected since small black holes have large 

r s than in the 
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Further, we have given several examples for the accretion flow problem, using different background metrics and 
equations of state, and analyzed cases where our assumptions are not satisfied. We have shown that violations of the 
assumptions on the fluid equation of state may lead to situations where the presence alone of a critical point is not 
sufficient to guarantee the existence of a globally defined accretion flow, while strongly distorted black holes violating 
our assumptions on the metric may lead to accretion flows which have distinct features compared to the standard 
case where all the assumptions are satisfied. 

There are various interesting questions that remain to be studied. For example, it would be desirable to consider 
more realistic fluids, where effects due to viscosity and emission of radiation are taken into account. Moreover, it 
should also be interesting to consider fluctuations away from the precise spherical and steady-state assumptions made 
in this work. Spherical and nonspherical linear perturbations of the Michel flow have been analyzed by Moncrief [7 
and shown to remain bounded outside the sonic sphere, thus establishing the linear stability of the Michel flow in the 
subsonic region. It would be interesting to extend this stability result to the region between the black hole horizon 
and the sonic sphere. In Ref. [52 we analyze small spherical and nonspherical acoustic perturbations of the Michel 
flow and show that these perturbations exhibit quasi-normal oscillations. 
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Appendix A: Justification for the metric conditions (Ml)—(M4) 

In this appendix we provide some additional details regarding the physical motivation for the metric conditions 
(M1)-(M4) given in Sec. [TTJ For definitions and discussions of Komar mass, Killing horizons, surface gravity and the 
energy conditions used below we refer the reader to Refs. [23] [25j. 

Condition (Ml) implies that the metric is asymptotically flat and that its Komar mass is well-defined and positive. 
The Komar mass is defined as 


MKomar — lim 
r—t oo 


87 tG 


-kdk. = lim 

r—>oo 


c 2 r 2 a'{r) 

2 G a(r) ’ 


where S r denotes a sphere of radius r, k = g(k, •) = g llu k^ 1 dx v = —acdt is the one-form associated with the Killing 
vector k, and * denotes the Hodge dual. According to our assumptions in (Ml) the limit exists and M^ omar = c 2 m/G. 
Next, condition (M2) implies that the Killing vector field k is null at the surface 77 := {r = rn}- Since a(r H ) > 0 
we also have N(th) = 0 which shows that the surface 77 is null. Therefore, 77 is a Killing horizon, and its associated 
surface gravity k, defined by da\ u = —2 c~ 2 n k| w is given by 

_ c 2 a'{r H ) 

K ~ 2 a(r H ) ’ 


which is strictly positive according to our assumption. 

Next, we provide justification for condition (M3), assuming that the stress-energy tensor T$) associated with the 
metric in Eq. ([Tj) via Einstein’s equations satisfies suitable energy conditions. The tt, rr and dd components of 
Einstein’s equations yield, respectively, 


1 



m 

N a’ 
r a 
TVcT 
r a 


47rG r 2 (0) 
c 4 ’ 

^ . 8 ^yo) 

r 3 ^ c 4 Vr 

8t tG (o) 

“u rPt , 


(Al) 

(A2) 


2or 2 dr 


(A3) 
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pi°^ and pj. 0 ^ refer to the energy 


where m(r) is the Misner-Sharp mass function, defined by N = 1 — 2 m/r, and 
density, radial and tangential pressure of the effective stress-energy tensor Tfp). (The superscript indicates that 
these quantities refer to the effective stress-energy tensor T^°J defined from the background metric as opposed to the 


stress-energy tensor T^ v of the accretion flow.) Combining Eqs. ( Al|A2|A3l we also obtain the following equations: 



N a’ 

4t tG 


r a 

c 4 

1 d . 


87 tG 

ar 2 dr 

c 4 


(^ (0) +^ 0) ). 
( e^+pW 


2 p 


(o) N 


(A4) 

(AS) 


The weak energy condition , stating the all physical observers should me asur e a nonnegative energy density, implies 
that > 0 and + p^ > 0. Therefore, under this condition, Eq. (Al) and the fact that m(ru) = th/ 2 > 0 
imply that the mass function m(r ) is positive for all r > rjj, and Eq. (A4) and a —> 1 for r — > oo yield a(r) > 0 


for all r > th- Under the additional assumption that the radial pressure is nonnegative, Eq. (A2) also implies that 
the function a is strictly monotonically increasing, justifying condition (M3). Instead of assuming pj 0 ^ > 0, it is also 
possible to resort to the strong energy condition , which states that e^ 


Pr 


^ + 2 p^ > 0. Eq. (|A5|) then implies 


that the function J(r) := r 2 a'(r)/a(r) is monotonically non-decreasing, which, together with J(rn) = 2r 2 H n/c 2 > 0 
implies that cr'(r) > 0 for all r > r#. 

Finally, we make some co mments regarding condition (M4). First, we remark that the strong and weak energy 
conditions and Eqs. ( A4|A5 ) imply that ( r 2 o')' > 0, such that the first inequality in Eq. |2]) is automatically satisfied. 
Next, consider the simple example in which a is a second-order polynomial in 1/r, 

, . 2 m em 2 

CT ( r ) = 1 _ _ + 

with e < 1 a parameter. The horizon is located at r# = m[ 1 + y/1 — e] and 

9 ,, . „ 2 em 2 

r"<7 (r) = 2m -> 0 


for all r >th- The condition (M4) is equivalent to 


T 9 

2e — 3— < 0 < --; 

m 8 + J 


f m\ 

1 - e— - e, r>r H , 
\ r 1 


which is satished for all negative values of e and all small enough positive values such that 

9 yJY^~e 


e < 


9 + 8\/l — e 


(1 + \/l — e). 


Appendix B: Proofs of Lemmata |l] [2] and [3] 

This appendix is devoted to the proofs of the technical results described in Lemmata [lj [2] and [3j 


1. Proof of Lemma [l] 

In order to prove LemmaJI] we first note that the conditions (M1)-(M3) imply that a is a well-defined, nonnegative 
function satisfying a(l) = 0 and lim^^oo a(x)/x = lim^oo <j(x)/(x 2 <j' (x)) = rn/^nr. Next, a short computation 
reveals that 




(. x 2 a’)’ 


dx \x 

Using condition (M4) we conclude from this that 
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, „ 9 d /a\ _ 9 m a 

1 + 3a > x 2 — I — ) >1-=-, 

dx \xJ 8 + — x 


(Bl) 
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where we have set m := m/rn- The first inequality is equivalent to xa' < 1 + 4a. We claim that, as a consequence 
of the second inequality, a(x) < x/(2 to) for all x > 1. If not, there must exist a point x* > 1 for which l/(2m) < 
a(x*)/x* < 2/(3 to) an d d(a(x*)/x*)/dx < 0 since a(x)/x —► 1/(2 m) as x —> oo. However, it follows from the second 
inequality in Eq. (|B1[) that 


dx \x 


© 


, 9ma(x*) 1 

> 1 - — > - > 0, 
8 x* ~ 4 


leading to a contradiction. Therefore, a(x) < x/(2m) for all x > 1. Using this result and again the second inequality 
in Eq. (B1) we find 


8a — 8a— 8a — 8a 1 

xa' > 1 +-> 1 + 


2a (4a+l) 2 


1 + 


8+2k 


16a + 1 


> 0 , 


and the lemma is proved. 


2. Proof of Lemma [2l 


(a) By condition (F3), w = <1/3. Integrating both sides from z to z\ with 0 < z < Z\ we obtain 


log 


v{z) 


1 

< - 
“ 3 


l°g (w) = lo § ( 


1/3 


1 / o 

from which the first inequality follows with 5 := v(z\)/z^ . For the second inequality, we use this result and 
the definition of the dimensionless sound speed v(z): 


9 log f(z) 


= v(z) 2 > 5 2 z 2/z , 0 < z < Zi. 


dlogz 

Integrating both sides from zq to z with 0 < zq < z < Z\ yields 


log 


m 

o) 


> -b 2 z 2 ^ 

- 2 


Taking the limit zq ► 0 and observing that f(zo) —> 1 by (FI), the second inequality follows. 

(b) From (a) we have f(z) > U/ 21 ) /3 for all 0 < z < z\. Taking the limit z\ —> z yields f(z) > > 1 

for all z > 0. Since f(z) — > 1 as z —>• 0, it follows that 


lim viz) = 0, 
z-v o 


as claimed. Next, it follows from (F3) that v is a monotonously increasing function which is bounded from above 
by one according to condition (F2). Therefore, when z —> oo, v converges to a finite value which is smaller than 
or equal to one. 

(c) By the definition of v and its monotonicity, we have for all z > zq > 0, 


^ log/ 

9 log z 


= v 2 (z) > v 2 (z 0 ). 


Integrating, we obtain from this 


f(z) > f(z o) 



z > ^o. 


where q := v(zq) 2 , which completes the proof of the lemma. 
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3. Proof of Lemma [31 


When 2 —>• 0, f(z) -A 1 and v(z) —> 0 according to assumption (FI) and Lemma |2](b), implying in particular that 
x c (z) —> 00 . As a consequence of condition (Ml) we obtain a(x c (zY) —> 1, and C(z) -> 1. On the other hand, when 
z->oowe have f(z) —> 00 and v{z) —> v\ < 1 according to Lemma[2](c) and (b). If iq < 1, x c (z) converges to a value 
larger than one, implying that cr(x c (z)) converges to a value larger than zero. Consequently, C{z) —> 00 as z —► 00 . If 
iq = 1 we have x c (z) —> 1 and a(x c {z)) —)• 0. However, in this case we can use L’Hopital’s rule and the relation 


to conclude that 
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dx c 

dz 
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1 w(z) 
2 a'(x c (z)) v(z) 2 


lim l = lim At 

2 ->oo 1 — v(z c y z->oo Aa'(x c (z)) zqz) 4 



(B2) 


Together with f(z) 00 this implies again that C(z) —> 00 . 

In order to pro ve th e nronotonicity statement we compute the derivative of C(z). Using the definitions of v and w 
and the relations (B2) and 


x c (z)a'{x c (z) _ 1 _ 4 v(z) 2 

cr(x c (z)) a(x c (z)) 1 -is(z) 2 ’ 


we obtain 
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1 - v(z) 2 


1 — v(z) 2 + w(z) — 


w(z) 

x c (z)a'(x c (z)) v(z) 2 J ’ 


2 > 0. 


(B3) 


Using Lemma [l] and 4 a(x c (z)) = v(z) 2 — 1 we find 


. , 1 w(z) , . 16a(x c (z)) +1 w(z) , , 0 , . . 

w{z) - _ ... ,y 2 > w(z) - . ^ = —3(1 - v(zf)w(z). 


x c (z)a'(x c (z)) v(z) 2 


[4 a(x c (z)) + l] 2 v(z) 2 


Since v(z) 2 < 1 and 0 < w < 1/3 it follows that the derivative of C(z) is strictly positive, implying its strict 
monotonicity. 
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